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1.  Introduction 


Long-lived  power  sources  offer  high  potential  in  applications  of  power  sources  for  unattended 
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sensors  and  communications  nodes.  ’  Recent  commercial  advances  have  been  described  with  the 
goal  of  developing  tritium  (  H)-powered  silicon  carbide  (SIC)  direct  energy  conversion  (DEC) 
sources  generating  energy  harvesting  levels  of  power.^’"^’^  In  a  direct  band  gap  material, 
absorption  of  light  is  generally  greater  because  there  are  fewer  competing  (losses)  mechanisms 
for  energy  transfer.  Common  direct  band  gap  materials  include  the  III-V  semiconductors  gallium 
nitride  (GaN),  gallium  arsenide  (GaAs),  and  indium  phosphide  (InP).  Common  indirect  band 
materials  are  silicon  (Si),  SIC,  diamond  (C),  and  germanium  (Ge).  Therefore,  direct  band  gap 
materials  can  be  more  efficient  at  elevating  valence  electrons  into  the  conduction  band  for 
photovoltaic  applications  than  indirect  band  gap  materials. 

A  series  of  simulation  were  designed  and  performed  to  evaluate  the  possibility  of  using  GaN  as  a 
radiation-tolerant,  wide  band  gap  semiconductor.  The  semiconductor  material  can  convert  the 
kinetic  energy  of  emitted  betas  in  radioactive  decay  emissions  into  a  series  of  lower  energy 
electron-hole  pairs  (EHPs)  and  knock-on  electrons.  These  secondary  electrons  are  swept  away  in 
the  electric  field  of  the  depletion  region  in  the  semiconductor  junction  and  collected,  forming 
current  flow  in  a  circuit.  The  nW/cm  power  generated  is  typical  of  energy  harvesting  levels  of 
power.''  Radioisotope  power  sources  differ  from  typical  renewable  energy/power  levels  in  that 
they  are  continuously  on  for  the  lifetime  of  the  isotope  decay.  The  power  output  degrades  by  half 
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with  each  half-life  (12.6  years  for  H)  of  the  isotope,  so  useful  energy  can  be  harvested  for  at 
least  two  half-lives  (decades). 


2.  Problem  Geometry 


The  simulations  model  a  volume  of  gas  emitting  a  beta  spectrum  from  ^H.  The  maximum  energy 
emitted  from  the  decay  of  is  18.6  keV.  The  decay  reaction  +p  +  v  creates  a  helium, 

beta,  and  anti-neutrino.  The  final  beta  spectrum  from  H  is  shown  in  Pig.  1.  12.5%  of  the  emitted 
betas  have  2.5  keV  energy.  The  average  P-energy  is  5.7  keV.  A  sample  GaN  p-i-n  device  was 
fabricated  by  metallo-organic  chemical  vapor  deposition  (MOCVD)  at  the  State  University  of 
New  York  (SUNY).  The  sample,  illustrated  in  Pig.  2,  was  modeled  using  the  Monte  Carlo  n- 
particle  extended  (MCNPX)  nuclear  scattering  code."*  MCNPX  is  a  general-purpose  Monte 
Carlo  code  that  can  be  used  to  model  neutron,  photon,  and  electron  (or  coupled)  transport.  The 
code  uses  an  extensive  collection  of  cross-sectional  data  and  is  able  to  simulate  the  transportation 
of  these  particles  with  energy  from  1  keV  to  100  MeV  in  materials. 
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Fig.  1  The  intensity  of  P-emission  spectrum  of  shows  that  the 
largest  P-flux  is  at  2.5  keV  (12.5%  of  the  flux) 


Fig.  2  Side  view  of  a  cylindrical  geometry  composed  of  GaN 
p-i-n  fabricated  for  initial  energy  conversion  evaluations 
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As  mentioned,  the  H  beta  emission  spectrum  is  shown  in  Fig.  1 .  ’  The  shape  of  the  beta 
spectrum  differs  in  isotopes  because  of  screening  potentials  and  forbidden  transitions.  Therefore, 
using  average  energy  or  endpoint  energy  is  not  sufficient  for  a  detailed  understanding  of  energy 
deposition  in  a  material.  The  shape  of  the  beta  emission  spectrum  is  responsible  for  variation  in 
energy  deposition  as  a  function  of  depth  in  materials.  The  H  energy  spectrum  is  defined  in  the 
source  definition  section  of  the  MCNPX  input  deck  detailed  in  the  Appendix. 

The  energy  emitted  in  a  decay  of  H  is  18.6  eV.  One  electron  is  emitted  in  each  decay  of  the  H 
isotope.  The  P  spectral  intensity  Ip  shows  that  12.5%  of  the  betas  are  ~2.5  keV.  The  median 
energy  is  6.1  keV.  The  P  spectrum  of  Fig.  1  is  used  as  input  for  the  MCNPX  calculation,  as 
described  in  the  source  definition  section  of  the  Appendix. 
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3.  Results 


The  electron  flux  is  calculated  in  MCNPX,  and  then  compared  to  a  model  of  EHP  creation 
efficiency  and  electron  beam  induced  current  (EBIC)  measurements  in  semiconductors.  The 
expected  values  of  number  of  EHPs  created  in  each  case  are  compared. 

3,1  Electron  Flux  in  GaN 

The  electron  flux/cm  was  calculated  in  200-nm  increments  throughout  the  depth  of  the  GaN 
cylindrical  volume,  which  extends  from  0-2.58  pm  along  the  axis  of  the  cylinder  (Eig.  3a).  The 
graph  of  electron  flux  throughout  the  depth  of  the  material  shows  charge  accumulation  along  the 
surfaces  (both  top  and  bottom).  The  P-spectrum  in  each  of  the  eight  200-nm  layers  is  shown  in 
Eig.  3b.  The  spectrum  in  the  top  (first)  200-nm  layer  shows  energy  content  out  to  18.6  keV  as 
expected  in  the  full  H  emission  spectrum.  The  P-spectrum  at  the  bottom  end  of  the  GaN  volume 
is  limited  to  1 1  keV  maximum  content.  The  flux  at  the  bottom  of  the  GaN  device  is  10”^  times 
reduced  (2.58  pm  thick). 


a)  b) 

Fig.  3  Flux  of  electrons  (per  incident  decay  beta)  generated  in  GaN  as  a  function  of  depth  in  the  material.  The 

flux  drops  by  a  factor  of  5000  front  to  back,  b)  Spectrum  of  electron  energy  as  a  function  of  depth  in  GaN 
is  similar  to  decay  spectrum.  The  maximum  electron  energy  decreases  from  18.6  keV  in  the  top  layer 
to  1 1  keV  in  the  eighth  layer  (1600  nm). 
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The  existing  H  power  sources  purchased  from  mb-microtec  Inc  provide  120  mCi/cm  from  the 
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rectangular  parallelepiped  H  -filled  vials.  The  P-fiux  resulting  is  4.4x10  p/cm  *s.  This 
corresponds  to  a  current  density  of  4.4  nA/cm  .  The  initial  nuclear  power,  P„uc,  from  decay  P 
emitted  (5.7  keV  Eavg)  is  then  4pW/cm^. 

3.2  Evaluation  of  Free  Carriers 

The  analysis  of  Klein  describes  the  energy  consumed  in  the  creation  of  an  EHP  (£ehp).  The 

o 

energy  required  to  generate  an  EHP  is  understood  to  include  three  loss  mechanisms.  They 
include  1)  band  gap  energy  (Eg),  2)  thermalization  loss  or  residual  kinetic  energy  of  electrons 
with  insufficient  energy  to  reach  conduction  band,  and  3)  phonon  losses  (lattice  vibrations,  r). 

(1) 


where  £ehp  is  the  energy  consumed  in  creation  of  an  EHP  and  r  represents  the  phonon  losses 
(Raman  quanta  emitted  by  ionization  impact).  The  9/5  factor  in  Eq.  1  is  experimentally 
determined  as  useful  to  describe  thermalization  losses. 


The  band  gap  of  GaN  (3.4  eV)  is  small  compared  to  the  Eavg  of  incident  P  decay  of  5700  eV.  The 
excess  energy  available  per  incident  P  could  create  a  maximum  of  (5700/3.4)  1676  secondary 
electrons  if  no  energy  was  lost  in  the  process.  However,  the  phonon  and  thermal  losses  reduce 
the  efficiency  of  secondary  creation  to  30%  in  GaN.  This  is  equivalent  to  requiring  -'lOeV  from 
each  incident  P  to  create  an  EHP  instead  of  just  the  band  gap  energy  of  3.4  eV.  Therefore,  the 
number  of  secondaries  created  could  be  as  many  as  570  per  each  incident  P  of  average  energy 
5700  eV.  Given  that  the  P  flux  calculated  for  100  mCi/cm^  ^H  vials  of  4.4x10^  p/cm^*s,  then  we 
might  expect  as  many  as  2.5x10  p/cm  *s  within  the  volume  of  the  intrinsic  region  of  the  GaN 
converter  or  400  nA/cm  .  How  well  we  can  collect  these  free  charges  from  the  device  is  another 
question. 


The  number  of  EHPs  generated  in  the  EBIC  diagnostic  technique  is 


Nehp  = 


3.2E 


ocE.  ^ 

■a — 


(2) 


where  E  is  the  energy  of  the  incident  electron.  Eg  the  energy  gap,  Ebs  the  mean  energy  of 
backscattered  electrons,  and  a  the  backscattering  coefficient.^  The  measured  values  for  number 
of  EHPs  created  in  the  EBIC  diagnostic  technique  has  been  experimentally  quantified  in  Eq.  2. 
Using  typical  values  for  backscatter  energy  of  10%  of  incoming  P  energy  and  a  backscatter 
coefficient  of  1,  we  can  calculate  an  upper  limit  for  the  number  of  EHPs  expected  as  a  function 
of  incoming  P  energy.  The  result  is  shown  in  the  top  curve  of  Pig.  4.  The  second  curve  in  Pig.  4 

is  calculated  using  the  weighted  number  spectrum  of  H  emission  from  100  mCi/cm  .  By  this 

12 

approach,  as  many  as  4.4x10  free  electrons  (area  under  curve)  could  be  available  for  collection, 

2 

possibly  generating  as  much  as  750  nA/cm  . 
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The  analysis  of  Bresse^  accounts  for  the  number  of  EHPs  created  from  incoming  electrons  while 

o 

the  analysis  of  Klein  focuses  on  the  mechanisms  of  loss  in  creation  of  EHPs.  Both  methods 
offer  an  upper  bound  of  free  carriers  in  the  device.  Experimental  measurement  of  EBIC 
techniques  can  generate  as  much  as  750  nA/cm  in  free  carriers  that  could  be  harvested  for 
potential  use  in  a  GaN  DEC  device.  By  calculating  the  energy  lost  in  the  creation  of  an 
individual  EHP,  as  much  as  400  nA/cm  could  be  generated.  These  expected  values  are  at  least  in 
the  same  order  of  magnitude. 


Energy!  keV)  keV 


Fig.  4  Number  of  EHPs  created  from  one  incident  electron  as  a  function  of  energy 
is  shown  compared  to  the  number  of  EHPs  calculated  from  100  mCi  of  the 
tritium  weighted  spectrum 

3,3  Electron  Energy  Deposition  in  GaN 

In  answer  to  the  question  of  how  big  should  one  design  the  i-layer  in  order  to  maximize 
maximum  charge  collection  efficiency,  the  MCNPX  simulation  numerical  result  of  energy 
deposition  in  the  material  provides  guidance.  The  MCNPX  simulation  models  nuclear  scattering, 
not  electron  transport  in  fields.  Therefore,  n-  and  p-doping  differentiation  is  not  significant  in 
these  calculations.  The  top  of  the  cylindrical  GaN  structure  is  uniformly  illuminated  over  the 
entire  surface  with  the  P  emission  spectrum  of  H,  as  shown  in  Pig.  2. 

The  logarithmic  contour  plot  of  energy  deposited  (MeV/g)  in  the  GaN  device  is  shown  in  Pig.  5a. 
The  mass  of  a  200-nm-thick,  1-cm-diameter  GaN  layer  is  96  pg.  Eight  200-nm  layers  are 
included  in  this  model.  The  two  dimensional  (2-D)  image  (Pig.  5a)  should  be  compared  to  the 
illustrated  line  drawing  of  the  device  in  Pig.  1.  The  energy  deposited  is  integrated  over  distance 
in  Pig.  5b.  We  see  that  68%  of  the  energy  is  deposited  in  150  nm  and  99%  of  the  energy  is 
deposited  in  700  nm. 
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a)  b) 


Fig.  5  a)  Cross-sectional  image  of  electron  energy  deposition  inside  the  GaN  DEC.  b)  Electron  energy 
deposition  is  integrated  along  the  depth.  68%  of  the  energy  is  deposited  in  150  nm. 


Corroborating  results  are  always  useful  to  track  when  using  large  simulation  codes  in  order  to 
confirm  that  the  correct  physics  models  were  chosen  for  the  code.  Gallium  characteristic  x-rays 
result  in  the  photon  energy  spectrum.  They  are  observed  at  the  K-alpha  and  K-beta  lines  of 
9,251.74  and  9,224.82,  respectively.^  These  low  energy  x-rays  do  not  contribute  in  a  significant 
way  to  any  radiation  dose,  as  only  81  of  lOOM  electrons  contribute  to  the  fluorescent  lines.  The 
simulation  also  predicts  that  0.5%  of  the  incident  P  contributes  to  a  bremsstrahlung  production 
spectrum,  also  an  insignificant  dose  level.  These  numbers  are  small  compared  to  24.2%  knock- 
on  electrons  generated.  The  current  generated  in  the  GaN  device  from  knock-ons  (electron- 
electron  scattering)  is  expected  to  be  1x10  eVcm  or  ~  160  pA. 


The  statistical  figure  of  merit  (FOM)  common  to  all  results  in  the  simulations  is  proportional  to 
the  square  of  the  signal-to-noise  ratio  (SNR).  The  FOM  is  recorded  for  each  signal  or  spectrum 
tally.  The  standard  deviation  of  each  numerically  calculated  result  is  determined. 


FOM  =  ^ 
<j^t 


(3) 


As  can  be  seen  from  Eq.  3,  the  FOM  is  inversely  proportional  to  t,  the  time  it  takes  to  run  the 
simulation,  which  is  directly  proportional  to  the  number  of  histories  calculated.  The  standard 
deviation,  a,  measures  the  uncertainty  related  to  statistical  fluctuations.  It  does  not  quantify  the 
accuracy  (or  truth)  as  it  relates  to  the  physical  quantity. 
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3.4  Comparison  of  Electron  Range 

The  MCNPX  calculation  described  in  Section  3.3  provides  more  detail  than  using  the 
continuously  slowing  down  approximation  (CSDA)  that  was  formalized  and  tabulated  into 
National  Institute  of  Standards  and  Technology  (NIST)  reference  tables.  The  three  energies  of 
significance  in  a  P-spectrum  that  describe  the  spectrum  shape  are  the  maximum  energy,  average 
energy,  and  energy  with  the  highest  flux  intensity.  For  H,  those  energies  are  18.6,  5.7,  and 

2.5  keV,  respectively.  The  CSDA  range  in  GaN  for  each  energy  is  shown  in  Table  1.  The  range 
at  which  68%  and  95%  of  the  free  carrier  energy  is  deposited  into  the  GaN  for  H  is  calculated  in 
the  MCNPX  simulation,  resulting  in  0.15  and  0.7  pm,  respectively.  It  is  interesting  that  the 
CSDA  range  (0.13  pm)  for  the  most  probable  electron  energy  (2.5  keV)  is  close  to  the  distance 
for  68%  energy  deposition  (0.15  pm)  and  the  CSDA  range  (0.7  pm)  for  average  energy 

(5.7  keV)  is  equal  to  the  range  of  95%  energy  deposition.  The  shape  of  the  beta  spectrum  will 
impact  this  correlation,  so  it  is  not  a  generally  useful  analogy;  however,  many  P  spectra  do  have 
a  similar  shape  to  the  (unlike  a  nickel  isotope  [^^Ni]  and  promethium  isotope  ['“^^Pm]). 


Table  1  CSDA  electron  range  in  GaN  for  three  commonly  available  isotopes 


NIST 

NIST 

NIST 

Eiviaxlnt 

eRange 

Eavg 

eRange 

Efnax 

eRange 

3H 

63Ni 

147Pm 

keV 

um 

keV 

um 

keV 

um 

2.5 

0.13 

5.7 

0.7 

18.6 

4 

1 

0.03 

17.2 

4 

67 

42 

1 

0.03 

62 

40 

224 

330 

4.  Conclusions 


H  is  the  most  available  and  inexpensive  of  the  isotopes  mentioned  in  this  report.  The  modeling 
result  using  the  full  P  spectrum  of  H  suggests  that  the  i-layer  should  be  between  150-700  nm 
thick  for  H  fuel  in  order  to  enclose  68-99%  of  the  generated  free  carriers.  The  optimized  ranges 
for  ^^Ni  and  Pm  would  be  larger  because  of  energies  typically  3  and  12  times  larger  than  that 
of  H.  The  optimized  depths  for  these  isotopes  will  be  calculated  in  a  future  report. 

In  these  calculations,  it  is  seen  that  the  use  of  energy  defined  by  maximum  flux,  the  average  over 
spectrum,  and  maximum  of  H  beta  spectrum  can  be  used  as  indicators  of  energy  distribution  in 
the  material;  however,  a  more  detailed  Monte  Carlo  analysis  provides  a  higher  degree  of 
accuracy  and  understanding. 
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Three  approaehes  to  estimating  the  upper  limit  of  the  free  earners  ereated  have  been  deseribed. 
The  energy  required  (and  loss  meehanisms)  in  EHP  ereation  is  largely  based  on  the  band  gap  of 
the  semieonduetor  and  phonon  losses,  and  results  in  a  400  nA/em  estimate.  The  number  of 
EHPs  generated  in  an  EBIC  diagnostie  teehnique  (based  on  Eg  and  energy  of  ineident)  suggests  a 
result  of  750  nA/em  .  The  number  of  knock-on  electrons  generated  in  the  collisions  of  incident  P 
with  GaN  atoms  is  estimated  to  be  -'200  pA.  The  MCNPX  simulation  experimentally  arrived 
probability  distribution  functions  and  theoretical  constructs  are  meant  for  10  keV  and  above,  but 
have  been  useful  down  to  1  keV.  Therefore,  when  using  MCNPX  to  calculate  free  carriers,  a 
value  of  1  keV  is  definitely  to  be  considered  a  lower  limit,  because  many  inelastic  scattering 
events  below  1  keV  occur  before  the  10  eV  events  begin  to  have  impact  on  the  generation  of 
secondary  free  carriers. 
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11 


The  input  deck  for  the  calculation  of  electron  transport  in  gallium  nitride  (GaN)  is  documented 
below.  The  file  name  is  GaNDECTr4.txt.  The  first  two  sections  connect  materials  to  volume 
geometries.  The  materials  for  each  volume  are  defined  in  the  third  section.  The  source  stimulus 
is  defined  in  the  fourth  section.  Type  of  physics  to  be  included  in  simulation  is  defined  in  the 
fifth  section.  Choices  of  results  displayed  and  tabulated  are  defined  in  sixth  and  seventh  sections. 
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Model  of  GaN-DEC  3H 


c - 

c  energy  deposition  modelling 


Q  + - 

- + 

C  1 

Cell  Cards 

1 

C  + - 

- -1- 

c  #  mat  density  surface  data 

3  4  -6.15 

-3 

imp:e,p  1  u=l 

SGaN  layer  1 

4  4  -6.15 

-4 

imp:e,p  1  u=l 

SGaN  layer  2 

5  4  -6.15 

-5 

imp:e,p  1  u=l 

SGaN  layer  3 

6  4  -6.15 

-6 

imp:e,p  1  u=l 

SGaN  layer  4 

7  4  -6.15 

-7 

imp:e,p  1  u=l 

SGaN  layer  5 

8  4  -6.15 

-8 

imp:e,p  1  u=l 

SGaN  layer  6 

9  4  -6.15 

-9 

imp:e,p  1  u=l 

SGaN  7 

10  4  -6.15 

-10 

imp:e,p  1  u=l 

SGaN  8 

11  9  -3.98 

-11 

imp:e,p  1  u=l 

SA1203  1st- layer 

12  1  -.00001 

- 

12  imp:e,p  1  u= 

=1  STr  vol 

c  air  encasing 


13  2  -0.00129  -13  345  6  7  8  9  10  11  12  imp:e, pi  u=l 

14  0  13  imp:e,p  1  u=l 

c  room 

15  2  -0.00129  -14  imp:e,p  1  fill=l 

16  0  14  imp:e,p0  Soutside  world 

c  + - -(■ 

c  I  Surface  Cards  | 

c  -1- - + 


c  #  type  param 

c  2  rcc  0  -.00001  0  0  0.00001  0  .5  $Si02  lOOnm 

3  rcc  0  0.0  0  0  0.00002  0  .5  $GaNp  200nm 

4  rcc  0  0.00002  0  0  0.00002  0  .5  $GaN  i 

5  rcc  0  0.00004  0  0  0.00002  0  .5  $GaN  n 

6  rcc  0  0.00006  0  0  0.00002  0  .5  $GaN  4 

7  rcc  0  0.00008  0  0  0.00002  0  .5  $GaN  n5 

8  rcc  0  0.00010  0  0  0.00002  0.5  $GaN  n6 

9  rcc  0  0.00012  0  0  0.00002  0  .5  $GaN  7 

10  rcc  0  0.00014  0  0  0.00002  0  .5  $GaN  8 

11  rcc  0  0.00016  0  0  0.00002  0.5  $A1203  3 
c  source 

12  rcc  0  -0.0002  0  0  0.0002  0  0.5  $Tr  (2um  thick) 
c  airto  encase  GaN 

13  rpp-1.0  1.0  -.001.001  -1.1  1.1  Sencasing  air 
c  room  dimensions 

14  rpp-1.2  1.2  -1.2  1.2  -1.2  1.2  Sroomair 


c  + - + 

c  I  Material  Cards  | 

c  -r - -I- 


c  m#  isotope  percent  &(newline)  Scomment 
ml  1000  1.0  $Tr  .OOOlg/cc 

c  ml  28063  1.0  $63Ni  9.552  g/cc 

m2  8016  .3  7014  .7  $air  0.00129  g/cc 

m4  31000  .5  7000  .5  SGaN  6.15  g/cc 

m8  14000  .33  8000  .67  $Si02  2.27  g/cc 

m9  13000.4  8000.6  $A1203  3.98  g/cc 


c  + - + 

c  I  Source  Definition  | 

c  -I- - -t- 


sdefpai=e  pos=  0  -0.00001  0  & 
erg=d4  axs=0  1  0  & 
rad=dl  ext=d3  & 
vec=0  1  0  dir=l 


sil  0  .5  Sradius  of  beam 
spl -21  1  Suniform  distribution 
si3  0  .00001  Slinear  distribution 
sp3  -2 1  0  Sprobabiiities 

si4  0  .001  .0025  .0038  .005  .0063  .0076  .01  .013  .015  .016  .019 
sp4  0  85  102  99  92  83  72  50  26  10  5  0  $3H 

Q  + - -I- 

c  I  Data  Cards  | 

Q  + - + 

mode  pe  Snhdtsa  these  not  needed 

PHYSiP 

PHYSiE 

nps  24000000  $  .lMh~.26min  90Mh~3hr 

c  -I- - -I- 

c  I  Tallies  | 

c  + - + 

fcl4  flux  averaged  over  the  cell  in  particles/cm^2 
fl4:e345  6  7  89  10 
el4  0  38i  .02 

fcl6  energy  deposited  averaged  over  cell  MeV/g 
fl6:p  3456789  10 
el6  0  38i  .02 

fc26  energy  deposited  averaged  over  cell  MeV/g 
f26:e  3456789  10 
e26  0  38i  .02 


c 

c  + - + 

c  I  Mesh  I 

c  - -I- 

c 


tmesh 

c  fiux=#/cm2  $dose=rem/hr  $pedep=MeV/cm3 
rmesh41  :e  flux  pedep  dose 
c  depth  profile 
cora41  -.6  58i  .6 
corb41  -.0002  98i  .0008 
corc41  -.6  .6 

endmd 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


2-D 

two  dimensional 

C 

diamond 

CSDA 

continuously  slowing  down  approximation 

DEC 

direct  energy  conversion 

EBIC 

electron  beam  induced  current 

EHPs 

electron-hole  pairs 

EOM 

figure  of  merit 

GaAs 

gallium  arsenide 

GaN 

gallium  nitride 

Ge 

germanium 

InP 

indium  phosphide 

MCNPX 

Monte  Carlo  n-particle  extended 

MOCVD 

metallo-organic  chemical  vapor  deposition 

Ni 

nickel 

NIST 

National  Institute  of  Standards  and  Technology 

Pm 

promethium 

Si 

silicon 

SiC 

silicon  carbide 

SNR 

signal-to-noise  ratio 

SUNY 

State  University  of  New  York 
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